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ABSTRACT 



We analyze the important formation processes for the Mg I 457.1 nm line. This line is an intercom- 
bination line and the source function is close to the local thermodynamic equilibrium (LTE) value. 
The strong coupling to the local temperature and the relatively high population of the lower level 
(the ground state of Mg I) makes this line an ideal candidate for temperature diagnostics in the lower 
chromosphere/temperature minimum region. Linking the temperature probed to an absolute phys- 
ical height is non trivial because of Non-LTE ionization. We analyze the Non-LTE effects and find 
that photo-ionization from the lower energy levels together with cascading coUisional recombination 
dominate the ionization balance. Taking properly into account the line-blanketing in the UV is essen- 
tial for obtaining the right photoionization rates. The identification of the main Non-LTE effects in 
the line allows us to construct a "quintessential" model atom, ideal for computationally demanding 
tasks, e.g. full 3D and/or time-dependent radiative transfer. Furthermore we analyze the diagnostic 
potential to temperature of this line in solar-like atmospheres, by synthesizing the line from a series of 
parametrized atmospheric models. These models have been constructed with fixed effective temper- 
ature, but with a variable heat term in the energy equation to obtain a chromospheric temperature 
rise at different heights. We conclude that the line has a significant potential in the diagnostics of the 
lower chromosphere temperature structure, especially for cooler atmospheres, such as sunspots. 
Subject headings: Sun: chromosphere - sunspots 



1. INTRODUCTION 

Deriving the solar atmospheric temperature structure 
from the emergent sp ectrum is a classical problem (e.g. 
Gingerich eral]|197lMVernazza et allllQSltlMaltbv et al.l 
19861: iFontenla et al.lll993[ ). It became clear quite early 
that the ID semi-empirical models would need an in- 
creasing temperature in the chromosphere, to explain the 
increasi ng radiation temperature with increasing opacity 
in lines (jAthav fc Thom£iill961h . The existence of such 
a static tempera ture rise in the real su n has, however, 
been questioned (jCarlsson fc Steinlll995D . Hence, it is of 
great importance to have a good temperature diagnostic 
for the lower solar chromosphere. 

One such candidate is the Mg I 457.1 nm line (from 
now on the 457.1 line). Note that we here use the 
common practice in astronomy to use the wavelength 
in air for naming lines, in vacuum the wavelength is 
457.2 nm. The 457.1 line arises from the forbidden 
transition 2p^ 3s^ ^ 2p^ 3s 3p , hence the source 
function should be tightly coupled to the local tempera- 
ture fe.g.lAthav fc House 196^:' Altrock fc Canfieldl fl97l 
iHeaslev fc Alie n 1980; Mauas et all 119881) due to the 
dominance of collisions in forbidden lines. The fact that 
the source function is highly coupled to the local tem- 
perature makes the 457.1 line ideal as a diagnostic of the 
temperature at the height of formation of the line. Early 
attempts to use t his line as a t empe r ature diagnostic in- 
clude t he work oflWhiteet a l. (1972^:'Alt rock fc Camion 



1972 ): 



19801) : 



Altrockl (11974): iRutten. (,1977): Heaslev fc Allen 
Mauas et al.1 (Il988l ). The most important result 



from this work is that the 457.1 line is formed at around 
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500 km height in typical quiet sun ID semi-empirical 
models, and hence the line will probe the temperature 
minimum region. Furthermore, the 457.1 line has also 
been used as a temperatu re diagnost ic for cooler models 
such as sunspot models (jLites et al.l 119871) . The mod- 
elling often gives a non-observed central reversal in the 
line core leading to the conclusion that the chromospheric 
temperature rise must be at a higher position. This re- 
sult has, however, been questioned (Mauas et al. 198^ 
due to possible uncertainties in the atomic model. More 
comprehensive Mg I mo del atoms have been p ublished 
like the 66-level atom of ICarlsson et al.l (|1992[ ) and the 
83 level atom of lZhao et al l (|1998D . The 457.1 hne has 
also been used to derive the temperat ure structure in so- 
lar fe atures such as m agnetic elements (iBriand fc Solankil 
1 19951) and flares (e.g. iMetcalf et al.lll990al lbD 

In this paper we present analysis and results of the 
temperature diagnostic properties of the 457.1 line. In 
12] we analyze the important formation processes of the 
457.1 line. As a result of this work we present a minimal 
realistic model atom (a quintessential model atom) for 
the 457.1 line in 5j3l In Sjljwe parametrize the model at- 
mospheres by solving the radiation hydrodynamics equa- 
tions with a prescribed heating per unit volume. The 
atmospheres are allowed to relax and the resulting at- 
mospheres are used to solve the Non-LTE problem for 
the 457.1 line. Finally we summarize our results in ij5] 

2. FORMATION PROCESSES FOR THE 457.1 LINE 

To understand the important formation processes for 
the 457.1 line we perform and analyse ID radiative trans- 
fer calculations. 

2.1. Atmospheric models 

For this analysis we use two serai-emp irical ID models, 
the FAL C model ()Fontenla et al.lll993l) and the Spot M 
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Fig. 1. — Termdiagram of the 66 level atom {left panel), the 27 level quintessential atom {center panel), and the 11 level test atom {right 
panel). To avoid confusion we have refrained from labeling levels with energies above 6.5 eV. 



model (jMaltbv et al.lll986D . These models represent an 
average quiet sun atmosphere and an average sunspot 
umbra atmosphere (for a sunspot with an umbral radius 
greater than 5" from the middle of the sunspot activity 
cycle), respectively. At the formation height of the 457.1 
line, magnesium is predominantly neutral in the sunspot 
model and mostly ionized in the FAL C model. This 
large difference in the ionization makes these two mod- 
els ideal for testing the formation processes in different 
atmospheric conditions. 



2.2. Model atom 

The atomic mod e l used in this analysis was published 
by iCarlsson et al.l (|1992l ). In short, this is a 66- level 
atom, including 65 Mg I levels and one Mg II level. The 
model atom include all levels up to n = 9 and also 
the 71 = 10 levels for the Mg I 2/ 3s 10s ^S, Mg I 
2/ 3s 10s ^S", and Mg I 2/ 3s lOp ^P°, since these lev- 
els have so big quantum defects that their energy is 
lower than the n = 9 energy levels with high 1 values. 
Furthermore, the model atom includes 65 bound-free 
transitions and 315 bound-bound transitions. Only one 
intercombination transition is included, 2p^ 3s^ <-> 
2p^ 3s 3p ^P°, which in fact is the transition giving rise 
to th e 457.1 line (note th at the data for this line is taken 
from i Mauas et al] ()1988D V Also note that there is no 
distinction between the singlet and the triplet systems 
for levels with I > 3, i.e. similar levels are considered 
to be internally Boltzmann populated. Finally, for some 
levels (most notably the upper level of the 457.1 line), 
the fine structure has been combined into "terms", with 
the individual levels assumed to be Boltzmann populated 
with respect to each other and the energy taken as the 
weighted mean energy with the statistical weight of each 
of the fine structure levels as weight factor. A termdia- 
gram of this model atom can be seen in Fig[TJ 

In these calculations we adapt a Mg abundan ce of 
7.53 on a logarithmic scale (|Asplund et al.l [2005[) . All 
other heavy element abundances are taken from the same 
source. 



2.3. Formation processes in the FAL C atmosphere 

Using the Mg model atom and the FAL C model 
atmosphere, we solve the equation of radiative trans- 
fer and statistical equilibrium using the MULTI code 
(Carlsson 1986). Background opacities from nearby lines 
are treated with a scattering term included in the source 
function in a two-level approximation. The line data for 
these lines were retr ieved from the Vienna Atomic Line 
Database (VALD) fPiskun ov et all 119951 : iKupka et all 
1999; Ryabchikova ct al. 1991). After solving the NLTE 
problem, we take into account the fine structure of the 
upper level of the 457.1 line by distributing the term pop- 
ulation over the sublevels solving the transfer equation 
again for the 457.1 transition. 

Since the bound-free edges of many of the important 
lower levels occur in the " UV line haze" , it is necessary to 
account for the resulting line blocking when calculating 
the bound-free opacity. This is done by using an opacity 
sampling method which reads a precomputed table of 
line opacity as fu nction of waveleng th, temperature and 
electron density (jCoUet et al.ll2OO50 . 

We test the NLTE effects in the line formation pro- 
cesses, by comparing the LTE solution with the NLTE 
solution. The emergent profiles in the two cases can be 
seen in the upper left panel of Figl2l The difference in 
the two line profiles is very small. The response-function 
to intensity given a c hange in temperature (from now on 
denoted Rai st, see iMagainI 119861 : iFossum fc Carlssonl 
[200l for details), is seen in the upper right panel in 
Figl^l There is not much difference between the LTE 
and the NLTE Rai,st, but this small difference is ac- 
tually a coincidence. Since the lower level in the 457.1 
line is the ground state of Mg I, the population density 
of this line is set by the ionization fraction of the atom. 
The fraction of Mg I atoms to the total number of Mg 
atoms is seen in the upper right panel of Fig[21 It is 
clear that the relative ionization in LTE compared to 
NLTE deviate 20% or more, but it turns out that the 
under-ionization at greater atmospheric heights to some 
extent cancel the over-ionization at lower heights in the 
atmosphere. We investigate these NLTE effect by com- 
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Fig. 2. — FAL C {upper panels) and Spot M {lower panels). Left panels: Emergent profiles in LTE {dashed) and NLTE {dotted). Right 
panels: temperature {black solid) intensity response functions to temperature calculated in LTE {blue dashed) and NLTE {blue dotted), and 
neutral fraction in LTE(red dashed) and NLTE {red dotted). The diamonds show where t^=1 in LTE and the triangles show similar results 
in the NLTE case. Note the difference in the Y-axis for the intensity, neutral fraction and response functions between the two atmospheres. 
Relative intensity is relative to FAL C continuum intensity in both cases. 

paring the rates into the Mg II level from each Mg I level, 
see FiglH As indicated by Figl3l the dominating effects 
in the ionization balance of magnesium is radiative ion- 
ization from the lower levels, (which dominate because 
of their relatively high population) especially the Mg I 
2p^ 3s 3p level (level number 3), which has signif- 
icantly larger photo-ionization crossection compared to 
other similar levels. Furthermore, the high energy triplet 
energy levels play an important role in the ionization bal- 
ance through coUisional recombination (jCarlsson et al.l 
[1992; Shimanskaya ct al. 2000). 

2.3.1. Testing the sensitivity of the ionization fraction 

To further test the dominating processes in the ion- 
ization balance of magnesium, we perform a series of 
test runs with different degree of realism in the radiative 
transfer. We test the two main NLTE effects, namely the 
photoionization and the coUisional recombination. For 
the first effect we run our model atoms with and without 
taking into account the "UV line haze" . For the second 
effect, we construct a model atom which does not con- 
tain the recombination ladder, essential for the cascading 
recombination. This is done simply by removing all en- 
ergy levels, except the Mg II level, with energy above 
6.5 eV, see Fig [T] right panel for details on the resulting 
11 level atom. The 11 level model atom is quite similar 
to the older models of the Mg I atom. The results of 
these tests for the FAL C and the Spot M atmospheres 
are seen in Fig|3]and can be summarized as follows: The 
largest effect comes from neglecting line blocking with 





2.0 


-'O 


1 1 1 1 1 1 1 1 1 








^ Radiative 


CTj 






+ Coll, triplet 


o 


1.5 




"o 






X Coll, singlet 




















1.0 






U 








d 








o 








a 








N 


0.5 






ioni 
















CP 


0.0 








-0.5 














20 40 60 








Level number 



Fig. 3. — Net ionization rates into Mg II from the different 
levels of Mg I at = 1 in the 457.1 line in the FAL C model. 
The symbol size is scaled to be proportional to the rate. Radiative 
rates {diamonds) , collisions from singlet levels {crosses) and triplet 
levels {pluses). 

a similar, but slightly smaller, change in the results if 
the cascading recombination is neglected. It is clear that 
the NLTE effects are significant, and it is very impor- 
tant that they are included in a realistic manner. If 
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Fig. 4. — neutral fraction in the FAL C {upper panel) and Spot M 
{lower panel) models with different level of detail included in the 
radiative transfer calculations. Note the difference in scaling on 
the Y-axis between the two panels. Line styles are defined in the 
figure. 

the NLTE effects are included, but with a lacking re- 
alism, the deviations from LTE are generally amplified. 
Finally, to understand the under-ionization high in the 
atmosphere, and the over-ionization in the lower atmo- 
sphere we turn our attention to the radiation field. If Ji, 
is bigger than Si, we have more photons absorbed than 
emitted and vice versa. Investigating the radiation field 
for the most important bound-free transition, we see that 
is larger than Si, where we have over-ionization and 
smaller than Si, higher in the atmosphere where we have 
under-ionization. 

2.4. Formation processes in the Spot M atmosphere 

Similar to the FAL C model we solve the equations of 
radiative transfer for the Spot M model. The resulting 
line profiles, LTE and NLTE, are seen in the lower left 
panel of Figl2l The difference between the two profiles 
is significant, especially compared to the similar FAL C 
profiles. The NLTE solution shows a pronounced in- 
crease in the central reversal in the line center of the 
line. This is basically caused by the higher geometric 
position of the monochromatic optical depth in NLTE, 
where the temperature is higher compared to the tem- 
perature at the height of formation in LTE. The Rai,st 
is much wider in NLTE compared to the LTE case. Fur- 
thermore, the Rai,st has contributions in NLTE from 
high regions where the temperature has a significant in- 
crease. The reason for the increasing formation height in 
the Spot M model is an increase in the fraction of neutral 
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Fig. 5. — Net ionization rates into Mg II from the different levels 
of Mg I at T,/ = 1 in the 457.1 line in the Spot M model atmosphere. 
The symbols are scaled to be proportional to the rate. Radiative 
rates {diamonds) , collisions from singlet levels {crosses) and triplet 
levels {pluses) are plotted. 
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Fig. 6. — Net ionization rates into Mg II from the different levels 
of Mg I at = 1 in the 457.1 line using the quintessential 27 
level mg atom in the FAL C model atmosphere. The symbols are 
scaled to be proportional to the rate. Radiative rates {diamonds), 
collisions from singlet levels (crosses) and triplet levels (pluses) are 
plotted. Compare to Fig|3] 

Mg, see Figl2]lower right panel. To understand this devi- 
ation in ionization fraction, we plot the net rates into the 
Mg II levels, see FiglH As for the FAL C model, we get 
big contributions from the lower levels of the atom, but 
also the higher excited levels contribute significantly to 
the ionization fraction, but with a negative contribution. 
This, together with the cascading collisional recombina- 
tion, gives much lower ionization fraction for the Mg I 
atom in the Spot M model. From Fig|4] lower panel. 
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Fig. 7. — Emergent profiles for the difli'erent model atmospheres are plotted {left panel). The temperature structure of these models are 
shown {right panel). Line styles are showing the different atmospheric models and the corresponding emergent profile. The atmospheres 
are denoted by the logarithm of the column mass (in cgs units, g cm~^) of the chromospheric temperature increase. 



one can see that the ionization fraction is less sensitive 
to the number of levels in the atom, but much more so 
for the treatment of line blocking. This is due to the 
reduced relative importance of the cascading recombina- 
tion, compared to radiative transitions, see FiglS] 

3. QUINTESSENTIAL MODEL ATOM 

From the above analysis in Sj2] we have identified the 
dominating processes in the formation of the 457.1 line. 
It is then a reasonable question if it is possible to con- 
struct a "quintessential" model atom (the smallest atom 
which reproduces the important physical processes) for 
this line. Since it is very important to get the ioniza- 
tion fraction correctly for this line, such a quintessential 
atom must contain a recombination ladder, which satis- 
factorily mimics the cascading recombination seen in the 
66 level model atom. 

Following the gene ral procedures outlined by 
iBard fc Carlsscml plOl . we construct a quintessential 
457.1 line model atom. First, we remove all singlet levels 
with energy levels above 6.5 eV. This is justified by 
the small contribution to the overall ionization fraction 
from these levels, as seen in Figl3] and in FiglH If two 
levels have a difference in effective quantum number 
Aq < 0.2 and at the same time are almost Boltzmann 
populated {Af3 < 0.1) in both in the FAL C model 
and in the Spot M model they are merged. The new 
model atom is tested in the two atmospheric models, 
and only accepted if the maximum relative deviation in 
the emergent intensity (compared to the full atom) is 
less than 25%. This procedure is repeated until all levels 
have been tested. The resulting model atom is reduced 
to a 27 level atom; see Figdl which is a reduction to less 
than half the size of the original atom. 

Further reduction is possible, but only with more lib- 
eral acceptance conditions, which would affect the real- 
ism of the model atom. The new model atom mimic the 
recombination ladder in the full model atom, by super 
levels as seen in Fig|6]and FiglSj 



The quintessential model atom was also tested against 
the origina l atom using a ser i es of dynamic atmospheres 
taken from lCarlsson fc SteinI (|1997l ). This test gives a de- 
viation between the lineprofiles of less than 3% between 
the two models atoms, and deviations in the bisector 
Doppler shifts of usually less than 10 ms~^ and always 
less than 30 ms~^. 

4. PARAMETRIZATION OF THE TEMPERATURE 
DIAGNOSTICS 

In the above analysis we have seen that the 457.1 line 
shows a strong central reversal in a cool sunspot atmo- 
sphere, but in a hotter quiet sun model there is no such 
reversal. The reason for the lack of reversal in the hot- 
ter model is the high ionization fraction in such models, 
which deplete the number densities of the lower level in 
the 457.1 line and put the formation height below the 
chromospheric temperature rise. From here on we focus 
on cool sunspot models, since the 457.1 line has a clear 
diagnostic potential for such models. We are especially 
interested in the possible bound the central reversal put 
on the position of the temperature rise in such models. 

To investigate the diagnostic properties of the line fur- 
ther, we construct a series of ID self consistent radia- 
tive hydrodynamics model atmospheres. We use the ra- 
dyn code (Carlsson & Stein 1992, 1995, 1997, 2002), to 
construct these model atmospheres. In short this code 
solves the mass, momentum, and energy conservation 
equations together with the detailed equations of radia- 
tive transfer for three of the important elements in the 
energy balance (H, C a, and He) on an adaptive grid 
(jPorfi fc Drurvl [l987f ). The boundary condition in the 
deep photospheric part (colum mass > 15 gcm~^) of the 
atmosphere is taken from a 4000 K effective temperature 
model without convection generated by the MARCS code 
(Gustafsson et al. 2008). The temperature structure in 
these deep parts of the atmosphere is fixed by adding a 
new term to energy equation in this region. 

The incoming UV radiation field from the corona has 
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to be specified at the upper boundary. T o facilita t e such 
a boundary condition we use the data oflTobiskal (ll99lD 
and follow the method used by |Wahlstrom fc Carlsson' 
(1994). Since the model used by Wahlstrom & Carlsson 
([1994) was based on a quiet sun corona, we expect that 
this model underestimate the coronal radiation field. We 
have performed tests (a more thorough analysis of this 
effect will be given in an upcoming paper) with different 
scaling of the incoming radiation field, and found that the 
radiation field has profound effects on the chromospheric 
temperature structure. From these tests we conclude 
that the incoming radiation field should be scaled up by 
at least a factor 5. With the above described boundary 
conditions, we allow the atmospheres to converge to an 
equilibrium state, and the pressure at the top of the sim- 
ulations is changed until there are no more outflowing 
material. The resulting temperature structures for an 
active and a quite corona are seen in FiglT] right panel. 
The difference in the chromospheric temperature struc- 
ture due to scaling of the coronal radiation field is quite 
big, however the effect on the 457.1 line is rather small 
due to the low formation height, see Figl7]left panel. For 
the construction of the atmospheres used in this analysis 
we adopt a scaling factor of 5 for the coronal radiation 
field. 

We parametrize the position of the chromospheric tem- 
perature increase in the theoretical model atmosphere by 
adding a constant heat term in the energy equation. To 
avoid unrealistic temperature structures in the lower part 
of the atmosphere, the heat term has to be defined per 
volume and not per mass. The latter will lead to an in- 
crease in temperature at very low atmospheric heights 
(|Bard fc Carlssonl[200l . Since the heating term per vol- 
ume will lead to a very high energy input in the corona, 
we allow a smooth transition to heat per mass above the 
transition region in the original atmosphere. The exact 
position of this smooth transition is changing the tem- 
perature structure in the corona, and to some extent the 
upper chromosphere, but it has no effect on the emer- 
gent 457.1 line profile. Furthermore, the heat term leads 
to an expanding corona but flows due to this expansion 
at the formation height of the 457.1 line are very small, 
typically less than 10 ms^^ and hence negligible. 

We construct different atmospheres by adding differ- 
ent heat terms per volume to the energy equation as 
described above. In this way we construct a grid of 
self consistent atmospheres with different position of the 
chromospheric temperature rise. Note that we use this 
parametrization of the temperature because it is only 
dependent on one parameter (the magnitude of the con- 
stant heat term added to the energy equation) . Further- 
more, these models should be considered a parametriza- 
tion of the position of the chromospheric temperature 
increase, and not a parametrization to evaluate the en- 
ergy balance in the chromosphere. 

Some of the resulting temperature structures are seen 
in FigOright panel. As seen in the left panel of Figl7]the 
emergent 457.1 line profile shows a significant increase in 
the central reversal with decreasing height of the chro- 
mospheric temperature rise. Note that these line profiles 
have been calculated using the quintessential modelatom 



and without including the nearby lines (to avoid confu- 
sion) . The size of the reversal is increased by the increas- 
ing temperature at the monochromatic optical depth of 
the line center of the 457.1 line. Note that the posi- 
tion of the monochromatic optical depth is only slightly 
changed with the changes in the position of the chro- 
mospheric temperature rise. The reason for this rather 
constant formation height is the high number densities of 
the Mg I ground state. If the temperature is increased, 
the relative change in the ionization fraction is changed 
only slightly due to the high density of Mg I atoms. 

Note the inconsistency between the emergent profile 
and the Eddington-Barbier approximation applied on the 
temperature structures in Fig[71 The reason for the lack- 
ing correspondence is the deviation of the source func- 
tion/temperature from a linear function in r^. Since the 
Eddington-Barbier approximation is not very good, it 
is not very instructive to invert the temperature pro- 
file from such assumptions. However, the reversals are 
strong enough to discriminate between different models. 
We recommend a forward modeling approach to the tem- 
peratu re diagnostics, similar to the work bv iLites et al.l 
(|1987t l. From these theoretical models the 457.1 line can 
be used to set lower bounds on the position of the chro- 
mospheric temperature increase. The sensitivity of such 
a discrimination will of course depend on the sensitivity 
of the observations used for the comparison. 



5. CONCLUSION 

We have analysed the dominant NLTE effects of the 
457.1 line in two solar models (FAL C and Spot M), and 
concluded that the ionization from the lower levels, espe- 
cially the Mg I 2p^ 3s 3p ^P°, is the main contribution to 
the ionization. Furthermore, cascading coUisional recom- 
bination in the triplet system plays an important role in 
the NLTE formation processes in the 457.1 line. With the 
understanding of the formation processes for the 457.1 
line, it is possible to construct a "quintessential" model 
atom. We have presented a 27 level quintessential model 
atom (an atom with the minimum number of levels, that 
maintain the important physical processes). 

We have found that increasing the realism in the model 
atom does not remove the non-observed central r eversal 
in the Spot M model, as claimed by Mau as et al.l ((1988). 

Finally we have presented a series of sunspot models 
where the height of the chromospheric temperature rise 
has been varied. These models show that the height of 
formation of the 457.1 line is not very sensitive to the 
height of the chromospheric temperature rise but the cen- 
tral reversal is. The amount of emission in the line core 
is thus a sensitive diagnostic of the lower chromosphere 
temperature structure and a lack of central reversal puts 
a strong constraint on the lowest possible position of the 
chromospheric temperature increase. 
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